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Abstract 

An instructive example is presented to elucidate the mathematical situation in the non¬ 
uniqueness problem of the infinite Friedmann-Keller hierarchy of equations for all multi¬ 
point moments within the theory of spatially unbounded Navier-Stokes turbulence. It is 
shown that the non-uniqueness problem of the Friedmann-Keller hierarchy emerges from 
the property that the system of equations is defined forward recursively. As a result, this 
system does not possess a unique general solution, even when the complete infinite system is 
formally considered. That is, even when imposing a sufficient number of initial conditions to 
this infinite system, it still does not provide a unique solution. This finding is supported by 
a Lie-group invariance analysis, in that the imposed example analogous to the Friedmann- 
Keller hierarchy admits an unclosed Lie algebra which allows for infinitely many functionally 
different equivalence transformations which all can be made compatible with any specifically 
chosen initial condition. Hence, if no prior modelling assumptions are made to close the 
Friedmann-Keller system of equations, the existence of an invariant solution within such a 
forward recursively defined system is then without value, since it just represents an arbitrary 
solution among infinitely many other, equally privileged invariant solutions. 

Keywords: Turbulence, Hopf Equation, Friedmann-Keller Hierarchy, Multi-point Correlations, 
Functional Equations, Ordinary and Partial Differential Equations, Infinite Systems, Lie Groups, 

Lie Algebra, Symmetries and Equivalences, Unclosed Systems, General Solutions, Cauchy Problems 
PACS: 47.10.-g, 47.27.-i, 05.20.-y, 02.20.Qs, 02.20.Sv, 02.20.Tw, 02.30.Mv 
MSC2010: 22E30, 34A05, 34A12, 34A25, 34A35, 40A05 

1. Introduction and motivation 

Due to the non-linear evolution of the Navier-Stokes equations, the statistical description of 
a turbulent flow on the level of the velocity moments is inherently unclosed. Presumably the 
first attempt to systematically formulate the equations of moments for turbulence was given 
by the proposed method of Keller & Friedmann (1924)F The result is an infinite hierarchy of 
coupled multi-point correlation (MPC) equations. Theoretical considerations show that on the 
level of the moments all (infinite) correlation orders have to be taken into account in order to 
allow for a consistent statistical description of turbulence, which, since then, became known 
as the still prevailing closure problem of turbulence. The motivation of Hopf (1952) was to 
formally bypass this problem in formulating a statistical description which operates on a single 
closed equation: Known as the functional Hopf equation^, it then induces the unclosed infinite 
hierarchy of MPC equations through a single moment generating functional <K 

‘Email address for correspondence: frewer.science@gmail.com 

tFor a more concise derivation of the Friedmann-Keller equations, see e.g. Monin & Yaglom (1971). 

*Note that also the discrete version of the Hopf equation, the Lundgren-Monin-Novikov hierarchy, which 
constitutes an infinite chain of coupled integro-differential equations, remains to be a closed system. The reason 
is that next to any possible initial conditions sufficient internal integral conditions are defined such that its 
general solution manifold is unique. For more details, see e.g. Hosokawa (2006) along with Appendix C in 
Frewer et al. (2014a) and the references therein. 
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Although one gained through this extension the theoretical advantage of dealing with a 
formally closed statistical description, the practical situation of constructing specific solutions of 
the problem, however, did not improve substantially when compared to the unclosed description 
of the moments. The fact is that on both levels of statistical description we face serious 
drawbacks: On the higher, more abstract level of the Hopf equation we face the problem 
of dealing with a functional calculus which from the outset is difficult to access numerically 
as well as analytically in a practical and constructive sense, while on the lower level of the 
Friedmann-Keller equations we face the problem of dealing with a hierarchy of coupled partial 
differential equations (PDEs) which is infinite. As I will show in this study by analyzing a 
lower-dimensional analogue, it becomes apparent that the latter problem faces a yet more 
fundamental problem: Even when formally considering all (infinite) equations in the statistical 
hierarchy, i.e. even if the Friedmann-Keller system is not truncated, the full system itself is still 
unclosed since no unique general solution can be constructed. In other words, when considering 
the case of a spatially unbounded flow e.g. in stating it mathematically as a Cauchy problem 
(as will be done throughout this study), the infinite system still is underdetermined and thus 
unclosed in that it does not offer any unique solutions although specific initial conditions have 
been posed. The reason for this defect is discussed in this study at the example of a lower¬ 
dimensional analogue, which shows all primary features of the Hopf-induced Friedmann-Keller 
hierarchy of MPC equations when stated as a functional Cauchy problem. These features are: 

i) the linearity of the higher level equation (the Hopf equation is a linear equation), 

ii) the linearity of the induced lower level moment equations (the infinite Friedmann-Keller 
hierarchy is a linear system), 

iii) the infinite hierarchy of moments is ordered forward recursively (each equation in the 
Friedmann-Keller hierarchy contains a higher order correlation function which only enters 
the next higher order equation in this chain). 

The non-uniqueness problem just mentioned above is rooted solely in feature iii), in that a 
certain hierarchy of differential equations is defined forward recursively. Because, as I will 
demonstrate, if the infinite chain of equations would be ordered oppositely, namely backward 
recursively (where each equation contains a lower order correlation function which then only 
enters the next lower order equation in this chain), then the non-uniqueness problem does not 
exist. But since the infinite Friedmann-Keller chain of moments is defined as a forward recursive 
system, it inherently faces the problem of non-uniqueness in its general solution manifold even 
if all infinite equations of the system are formally taken along. Consequently, when performing 
an invariance analysis on such a system, one can generate infinitely many and functionally 
diverse Lie-point symmetries^, and hence infinitely many different and independent invariant 
solutions which all can be made compatible to any posed initial condition. 

My motivation for this study was to elucidate this non-uniqueness problem for the MPC 
equations in all its facets, because it seems that in the relevant literature on turbulence 
there still exists a misconception on this issue, in particular in the studies of Oberlack et al., 
e.g. in Oberlack & Rosteck (2010), Oberlack & Zieleniewicz (2013), Avsarkisov et al. (2014), 
Waclawczyk et al. (2014), and Oberlack et al. (2014). In all these studies the reduced (lower 
level) Friedmann-Keller hierarchy of multi-point moments is incorrectly treated as a closed 

'In order to warrant consistency within an invariance analysis as for the Friedmann-Keller system, the actual 
symmetries must be identified as weaker equivalences. The reason is that since the considered system of equations 
is unclosed, all admitted invariances only act as equivalence transformations and not as symmetry trans¬ 
formations. For more details, see e.g. Section 2 in Frewer et al. (2014a), Sections 3-5 in Frewer et al. (20146), 
and the references therein. 
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system, with the misleading argument that it’s due to that the system is infinite dimensional, 
because ”... if the infinite set of correlation equations is considered the closure problem is some¬ 
what bypassed” (Oberlack & Rosteck, 2010, p. 454). Also the performed invariance analysis in 
each of these studies is misleading. For example, to obtain the result that ”... the latter set of 
symmetries is considerably enlarged for the infinite set of MPC equations” (Oberlack & Ziele- 
niewicz, 2013, p. 12)1 j s not a surprising result, because the considered infinite Friedmann-Keller 
system of moments, as we will assert herein due to its forward recursively organized hierarchy, is 
simply an unclosed and thus an underdetermined system of equations which intrinsically allows 
for arbitrary invariances. Furthermore, to address the open question that ”... so far complete¬ 
ness of all admitted symmetries of the MPC equation has not been shown... [which] appears to 
be necessary not only from a theoretical point of view but rather essential to generate scaling for 
all higher moments” (Oberlack et al., 2014, p. 1702), and that ”... hence, finding new statistical 
symmetries and/or proving the completeness of the set of symmetries is a next task for further 
study” (Waclawczyk et al ., 2014, p. 10), would collectively only lead to a superfluous task, 
because, as already said, arbitrary invariances can be constructed when considering an un¬ 
closed infinite hierarchy of equations as the MPC equations considered by Oberlack et al.; 
therefore completeness in this very sense can never be obtained in such infinite systems. 
Finally, when trying to correctly interpret any results arising from such an invariance anal¬ 
ysis, in particular when trying to construct invariant solutions for such unclosed systems, it 
is necessary to recognize the subtle but important difference between a symmetry transfor¬ 
mation and an equivalence transformation. While a symmetry transformation always maps a 
solution to another solution of the same (closed) equation, an equivalence transformation, in 
contrast, generally only maps a possible solution of one underdetermined (unclosed) equation 
to a possible solution of another underdetermined (unclosed) equation (see e.g. Ovsiannikov 
(1982), Meleshko (1996), Ibragimov (2004), Bila (2011)). Hence, since any Lie-group invariance 
analysis on the infinite Friedmann-Keller hierarchy will be based on equivalence and not on 
symmetry groups, it is misleading and even ill-defined to construct invariant solutions of the 
unclosed Friedmann-Keller equations if no modelling assumptions on these set of equations are 
being priorly invoked (Frewer et al., 2014a, b). Therefore, to state that ”... a variety of classical 
and new scaling laws” were derived for which ”it was shown that they are exact solutions of sym¬ 
metry invariant type of the infinite dimensional series of MPC equations” (Avsarkisov et al., 
2014, p. 102) is more than misleading, because what to understand under ’’exact solution” 
when the construction process for generating solutions itself is arbitrary? Independent of 
whether a Lie group invariance analysis is employed or not, these ’’exact solutions” are nothing 
else but arbitrary and thus non-privileged solutions. They just arise from an unclosed infinite 
dimensional system which does not give any indication of whether these ’’exact results” are 
relevant solutions or not. For that, additional, external information is needed, e.g. as posing 
certain modelling assumptions to close the infinite system of equations. 

The following study is organized as follows: Section 2 introduces the functional Hopf equa¬ 
tion in physical space, from which the infinite Friedmann-Keller system of multi-point moments 
for spatially unconfined turbulence is then reduced. Section 3 investigates a lower-dimensional 
analogue to the system established in Section 2, but, which first will only feature the proper¬ 
ties i) and ii) listed above. Instead of property iii), a backward recursively infinite system is 
first considered, in order to demonstrate that for such a case the non-uniqueness problem as 
discussed above cannot arise. Section 4 forms the key part of this study, in that an extended 
lower-dimensional analogue is considered which now will feature all properties i)-iii) of the 

'In Frewer et al. (2014a, b) it has been shown that this newly enlarged set of symmetries is even unphysical in 
that it leads to various inconsistencies with the underlying deterministic Navier-Stokes theory. Hence, besides 
the well-known classical symmetries of the Euler and Navier-Stokes equations (Fushchich et al., 1993; Olver, 
1993; Frisch, 1995; Andreev et al., 1998), no new enlarged physical set of statistical symmetries exist. 
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system established in Section 2. This analogue will allow us to elucidate step by step all facets 
of the non-uniqueness problem which inherently characterizes the infinite Friedmann-Keller 
system of multi-point moments. Section 5 concludes and completes the investigation. 


2. The Hopf equation and its induced MPC equations 


The functional Hopf equation (Hopf, 1952; McComb, 1990) for the characteristic or moment 
generating functional <f> = <J>[y(x), t] in physical space is given by the linear differential equation 


<9$ 

~di 


= J Va ( x ) M 


d 


5 2 <5> 


dxp Sy a (x)Sy^(x :) 


+ v A 


5 4> 


5y a (x) 


d 3 x, 


( 2 . 1 ) 


where y(x) is a real, integrable and time-independent auxiliary vector held being comple¬ 
mentary to the ensemble of all admissible incompressible velocity fields u(x) governed by the 
Navier-Stokes equations sampled at each time step t > 0. Since we consider an unbounded 
how, we assume that the ensemble of helds y(x) is vanishing sufficiently fast at infinity. The 
helds y 1 - in equation (2.1) represent the transverse (solenoidal) part of y(x) in order to elimi¬ 
nate the pressure terms which otherwise would arise in (2.1). We will use the transverse form 
as used in Hopf (1952) (see also Appendix A) 

y«(x) = y a (x) - d a cp(x), (2.2) 


with the scalar held 


(p(x) 


V ' • y( x 0 ,y 

47r|x — x'| 


(2.3) 


Note that in order to warrant physical consistency, a mathematical solution of the Hopf equation 
(2.1) is only admitted if for all times the following conditions for the characteristic functional 
are fulfilled 


4>*[y(x),t] =4>[-y(x),t], 4>[0,f] = l, |<h[y(x),f]| < 1. (2.4) 


Now, when reducing the MPC equations from (2.1), one can either work in the representation 
of the full vector held y, or in the decomposed representation of its transverse part y^. The 
former representation of the Hopf equation (2.1) takes the form (see Appendix B.l) 


<9$ f . J. 8 5 2 <S> . 4<I> \ )3 

dt J Va X \ dxp 5y a (x)5yi3{x) +U 5y a (x) J 

r / <9 1 \ d 2 _ 

J ^ a \dx Q 4 tt\x — x'\ J dx'^dxiy 8y/s(x')6y 1 (x') 

while the latter one takes the more concise form (see Appendix B.2) 


d 3 x' d 3 x, 


(2.5) 


<9$ 

~di 



(. d 5 2 <f> \ 3 

V dx p 8y^(x)5yj(x) +l> 8y^(x)) 


( 2 . 6 ) 


Note that both equations, (2.5) and (2.6), refer to the same characteristic functional 4*, since 
under the dehning relation (2.2) it transforms invariantly 4>[y(x),f] = < f , [y _L (x), t] (Hopf, 1952, 
p. 93). Keep in mind that equation (2.6) is a projection of equation (2.5) onto solenoidal vector 
helds, where if V • y ^ 0 then V • y^ = 0. In other words, if the ensemble of ah helds y would 
be divergence-free, i.e. V • y = 0 (which in general is not assumed), then (2.5) is identical to 
(2.6), as can be readily observed in (2.5) after executing a partial integration with respect to x. 

Further note that both representations (2.5) and (2.6) of the Hopf equation (2.1) lead to 
the same set of MPC equations. It’s obvious that the latter representation (2.6) is suited best 
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for any formal investigations as we are interested herein. It induces considerably simplified 
formal expressions for the MPC equations since no explicit pressure terms would surface as 
in representation (2.5). The corresponding pressure-dependent MPC equations, when needed, 
can then either be derived by using the defining relation (2.2) in the obtained results from 
the transverse-projected equation (2.6), or by directly deriving these from the full-field equa¬ 
tion (2.5). 

According to (Hopf, 1952, p. 100-101) the infinite Friedmann-Keller hierarchy of MPC 
equations can be generated from (2.6) by performing a Taylor expansion of t] around the 

functional point y^ = 0 : 


M yM 


y x =o 


+ 


y x =o 


U[y\(] = l + /rfW,( Xl )^ i(xi) 

+ ^ Kt (X!) 

which formally can be written as 

$[y\t] = l + $ 1 [y ± ,t]+$ 2 [y ± ,t] + ... , 

where ^"[y^t] is a homogeneous polynomial functional of degree n in y^ = y J_ (x), 
= j d 3 xi---d 3 x n2 /^ 1 (xi)--- 2 /^(x n )/\^ 1 Qri (xi,...,x n ,t), 
and where the components of the kernel function = (K^ a ) are given by 


^ai...a„( x l »••• >**».*) = -7 


(^[yV] 


n\ Sy^ (xi) • • • Sy^ n (x n ) 

5^[y,t] 


y x =o 


n! <5y Q1 (xi)---5y Q „(x n ) 


(2.7) 


( 2 . 8 ) 


(2.9) 


— A-aa(xi,..., x n , f). (2.10) 


y=o 


The last equality stems from the fact the functional derivatives stay invariant under the trans¬ 
formation (2.2) (see Appendix B.2); and since y^ = 0 o y = 0, the kernel function constitutes 
itself as an invariant function, i.e. = K n . In Hopf (1952) this function is determined as 
the equal-time MPC function of the incompressible Navier-Stokes velocity field u(x, t) 


i n 

K ai ...a n (x 1 , . • • ,Xn,t) = — (« ai (xi,i) • ■■U an (x n ,t)), (2.11) 

where (•) denotes the statistical ensemble average of the instantaneous and thus fluctuating 
velocity field evaluated at n different points. Upon inserting the power series of T (2.8) into 
equation (2.6) and upon equating terms of equal degree on both sides, one finally obtains the 
infinite sequence of differential equations for the MPC functions 


<9T n 

~df 


Va( x ) 


d 5 2 $ n+1 
l dx/3 <%K x )<fys ( x ) 


+ nA 


5$ r 


^Va ( x ) 


d 3 x, Vn > 1. 


( 2 . 12 ) 


By taking the notation from Oberlack & Rosteck (2010), the above system can be evaluated to 
(see Appendix C.l) 


mi 

dt 


+ E V *. 


H 


z,n+1 


-"Y, A 


x,H £ = 0 , Vn > 1, 


(2.13) 


i =1 


i— 1 
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where 

H^(x 1; ... ,x n ,t) = H n ( Xl) .. .,x n ,t) = (u(xi ,t) <8> • • • <8> u (x,t)), (2.14) 

is the instantaneous (equal-time) multi-point velocity correlation function of n-th order, and 
where 

Hy +1 (xi,... ,Xj,... x n , Xi,t) = ( lim = H n+ i(xi,... ,x n ,x n+1 ,t) ) , 1< 

’ \x n +l —j 

is not a (n + l)-point function, but only a lower dimensional n-point function 
order with the full-ranked transverse (solenoidal) property (see Appendix C.2) 

(V X1 (8) ■ ■ ■ <8> V Xi (8) ■ ■ ■ ® V Xn <8> V Xi ) • H^ n+1 (xi,... ,Xj,.. .x n ,Xj,f) = 0, 1 < 

The MPC equations (2.13) go along with the following natural incompressibility 

V Xfc • H n (xi ,... ,x n ,t) = 0, \/k between 1 < k < n, 

(xi,..., x;,... x n ,xi,t ) = 0, \/k, l between 1 < (, k , l ) < n, for k ^ 

due to the incompressibility constraint V Xfc • u (xf.,t) = 0 of the instantaneous 
u(xfc,i), when evaluated at each point x = x^ (for all k = 1,n) within the single physical 
domain x. Note that, when decomposing the transverse fields (2.15) into their full fields (see 
Appendix C.2), the system (2.13)-(2.17) matches the infinite Friedmann-Keller hierarchy of 
MPC equations as derived in Oberlack & Rosteck (2010) for the instantaneous case. 

Now, if we collect the key properties of the combined system (2.1) and (2.13), the most 
important thing about it is that it is linear. Both the higher level Hopf equation (2.1) as well as 
its induced lower level system of MPC equations (2.13) constitute a linear system of equations. 
The former as a singly closed functional equation, while the latter as an infinite hierarchy of 
PDEs. As we will see during this study, a crucial aspect of the considered hierarchy (2.13) 
is that it is defined forward recursively. To illustrate the basic problems from which such an 
infinite system suffers, we will investigate the following lower-dimensional analogous system for 
a spatially one-dimensional field u = u(x,t ) with its induced moments u n = u n (t ) 

dtu = d^u — A • x 2 u, (2-18) 

du f °° 

—2- = n ■ (n — 1) • u n -2 — A • u n+ 2 , where u n (t) = x n ■ u(x, t ) dx, n > 0, (2-19) 

(It .7 — 00 

which features all key properties of the original combined system: The linear and non-autono- 
mous property of the higher level functional Hopf equation (2.1) is featured by the higher level 
PDE (2.18), from which the lower level ordinary differential (ODE) moment equations (2.19) are 
induced, which themselves again feature the linear, autonomous and forward recursive property 
of the infinite Friedmann-Keller hierarchy of MPC equations (2.13) which again emerge from 
the higher level Hopf equation^ (see Appendix D for the derivation of (2.19) from (2.18)). 

When considering in particular the original system of moments (2.13), the foremost aim, 
of course, is to find certain solutions, at least some asymptotic solutions, e.g. as attempted 
by Oberlack et al. in applying a Lie group based symmetry group analysis directly on system 
(2.13) in order to generate certain specific invariant solutions which then should function as first 

^Note that the incompressibility constraints (2.17) cannot be mapped down to a one-dimensional system. 
Hence, the incompressibility property of the original MPC system of equations (2.13) is not featured by its 
analogous infinite system (2.19). But this is no drawback for the issues to be addressed in this study, since the 
incompressibility constraints are irrelevant and thus do not influence the conclusions made herein. 


i < n, (2.15) 
of (n + l)-th 

i < n. (2.16) 


constraints 



velocity field 
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principle scaling laws within some asymptotic flow regime (Oberlack & Rosteck, 2010; Oberlack 
& Zieleniewicz, 2013; Avsarkisov et al., 2014; Waclawczyk et ai, 2014; Oberlack et al ., 2014). 

Our claim is that system (2.13), along with its incompressibility constraints (2.17), is un¬ 
closed, even if all infinite equations are formally considered. Hence such an approach as per¬ 
formed by Oberlack et al. is misleading when a symmetry analysis is only carried out for the 
statistical system of moments without additionally incorporating the underlying deterministic 
set of equations. In other words, to generate any kind of solutions directly from (2.13) is not 
well defined if no prior modelling assumptions in accord with the deterministic Navier-Stokes 
equations are made to close the statistical system of moments. 

The reason why system (2.13) is unclosed, even if all infinite equations are considered, is 
basically twofold: The first issue lies within the non-symmetric and thus incomplete coupling 
between the equations. For each order, the unknown moment H- L ri+1 is not directly coupled 
to the next higher order equation. Although all components of that moment can be uniquely 
constructed from the higher dimensional moment H n+ i once they are known, which is formally 
denoted as (2.15), the necessary inverse construction, however, fails. Hence, since (2.15) is a 
non-invertible construction, i.e. since H n+ i cannot be uniquely constructed from H^ n+1 , these 
latter moments are to be identified as unclosed functions in system (2.13) as they do not directly 
enter the next higher order correlation equation. For each order n in the hierarchy the total 
number of dynamical equations (2.13), along with the continuity constraint equations (2.17), is 
therefore always less than the total number of unknown functions. In total this just reflects the 
classical closure problem of turbulence for the moments which cannot be bypassed by simply 
establishing the formal connection (2.15) — for a more detailed discussion on this issue, please 
refer to Appendix A & C in Frewer et al. (2014a) and Section 3 in Frewer et al. (20146). 

The second and more fundamental issue of why the infinite hierarchy (2.13) is unclosed, 
and which forms the focus of this study, is that the hierarchy is defined forward recursively. 
As I will elucidate at the example of the lower-dimensional analogue (2.19) to system (2.13), 
an infinite forward-recursively defined hierarchy of equations always leads to a non-uniqueness 
problem in its solution manifold — even when sufficient initial conditions are posed, a system 
such as (2.13), when stated as a well-posed initial value problem, still does not return a unique 
solution. However, this problem does not arise if any infinite hierarchy is considered which 
is defined backward recursively. Hence, instead of the full system (2.18)-(2.19), I will first 
investigate the solution properties of the subsystem A = 0, which will be done in the next 
section, thus leading us to an infinite hierarchy of moment equations which first will be defined 
backward recursively. 


3. Example for an infinite backward differential recurrence relation 
— A closed system with a unique general solution manifold 


In this section we will analyze the combined system (2.18)-(2.19) for A = 0, which hence results 
to a study of the usual one-dimensional diffusion equations and its moments. Consider first 
the higher level PDE (2.18) by stating it as the following well-posed Cauchy problem 

dtu = d^u, with u(x,0) = (3-1) 

where we assume that the function u = u{x , t) is decaying sufficiently fast at infinity, i.e. 


lim u(x, t) = 0, lim d x u(x, t) = 0, Vi > 0, 

>-d=oo x —>-=boo 


(3.2) 


and where, for convenience, we want to assume that the initial function (f is normalized to 

/ OO 

4>(x) dx = 1. 

-oo 


(3.3) 


Then the initial value problem (3.1) has the unique solution (see e.g. Polyanin (2002)) 


l roo ( x _ x ')2 

u(x,t) = —== / e « <f>(x')dx', for t > 0. (3.4) 

V47Tf J -oo 


If we would choose as a initial condition, for example, the Gaussian distribution 


4>(x) = — i=e 2 ( ) , cr, fi € M, u > 0, (3-5) 

cry 2vr 


then the solution (3.4) would read 


1 1 (:e-M) 2 

u(x,f) = ==e t > 0, (3.6) 

V ' V2vr(2t + (r 2 ) V ' 


with its associated moments uniquely given as (shown here only up to third order) 

POO 

n = 0 : uo(t) = x° ■ u(x, t ) dx = 1, 

J — OO 

(3.7) 

n = 1 : 

Ul{t) = j 

POO 

' x 1 ■ u(x, t ) dx = fi, 

— OO 

(3.8) 

n = 2 : 

u 2 (t) = j 

POO 

1 x 2 ■ u(x,t) dx = 2t + /r 2 + a 2 , 

—OO 

(3.9) 

n = 3 : 

«3(t) = J 

f x 3 • u(x, t)dx = n ^6t + fi 2 + 3a . 

— OO ^ 

(3.10) 


Now, let’s consider the corresponding moment equations when they are independently derived 
as an infinitely coupled hierarchy from the higher-level Cauchy problem (3.1) (see Appendix D) 

du poo 

— -p- = n ■ (n — 1) • u n ~ 2 , with u n ( 0) = / x" • <f>(x) dx, n > 0. (3-11) 

dt J —oo 


Our aim is to investigate in how far the above infinite system of first order ODEs represents a 
well-posed initial value-problem for the moments u n = u n (t). The first positive observation is 
that to each unknown function u n one can bijectively associate an initial condition u n ( 0) to it, 
since any first order ODE only requires a single initial condition in order to provide a unique 
solution. The next step is to derive the general solution, for which it is helpful to recognize 
that since the infinite hierarchy in (3.11) is defined backward recursively of an order two, it can 
be rewritten into the following equivalent form 


du 0 (t) 

dt 

du\{t) 


dt+ dt*+l 2(t) = ( 2t + 2 > ! ' “»(*)' *>0, 1 
^H+lW = (2i + !)!•„,(*), ;>1. J 


(3.12) 


dt 
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This system can then be uniquely integrated to give the general solution 


uo(t) 
U2k+2(t ) 


= c 0 , 


k (!) 

C2k+2 + E ~T tJ 
j =1 


rtk=t rtk -1 Ho 

+ {2k + 2)\ / / ••• / uo(t) dt'dto ■ ■ ■ dtk-i, k>0, 

Jo Jo Jo 


ui(t) = Cl, 


> (3.13) 


l-l ( 2 ) 

U2l+l(t) = C21+1 + E ~^f~ ^ 
j = 1 

rti=t rti _! Hi 

+ (21 + 1)! / / • • • / mi(£ / ) d£i • • • dtq_ i, l > 1, 

Jo Jo Jo 


with the expansion coefficients given as 


o (1) 



(2k+ 2)1 

(2k + 2 — 2j)\ C2k+2 ~ 2 i 


k > 0; 0 < j < fc, 


(21 + 1 )! 

(21 + 1 - 2j)\ C2l+1 ~ 2 i 


l > 1; 0 < j < l - 1, 


\ 




(3.14) 


where all c n for n > 0 are arbitrary integration constants. Hence, the second positive observa¬ 
tion is that the unrestricted system of equations in (3.11) 

= n ■ (n - 1) • u n - 2 (t), n > 0, (3.15) 

provides a general solution (3.13) which only involves arbitrary constants, i.e. the unrestricted 
system (3.15) provides a unique general solution and therefore it ultimately represents a for¬ 
mally closed system of equations. Because, when restricting this system to the underlying 
PDE’s initial condition u(x, 0) = <$>(x) as given in (3.1), which for the infinite ODE system 
takes the integral form u n ( 0) as given in (3.11), will turn the general solution (3.13) into a 
unique and fully determined solution, where the integration constants are then given by 

c n = u n ( 0), n > 0. (3.16) 

By choosing for example again the specific initial condition (3.5), the general solution (3.13) 
will give exactly the expressions (3.7)-(3.10) for the moments, which, as we know, were directly 
derived from the higher-level PDE solution (3.4) of the well-posed Cauchy problem (3.1). 
Hence, since the unrestricted infinite system of ODEs (3.15) constitutes a formally closed 
system of equations, the restricted system (3.11) expresses a well-posed initial value problem. 

That (3.15) constitutes a formally closed system is also supported when performing a Lie- 
group based symmetry analysis (see e.g. Stephani (1989); Olver (1993); Bluman & Kurnei 
(1996); Ibragimov (1994); Hydon (2000)). The infinite and coupled system (3.15) admits the 
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following complete and unique set of Lie-point symmetries^ 


Si: X 1 = 8 t , 

s 2 : X -2 = tdt - 2 u 0 d uo - u 2 d U2 -|-b (n - 2 ) ■ u 2 n d U2n H- 

2u\d ui u^d U3 + • • • + (n 2) • it 2 n-bi^?i 2 n+i T ''' > 

5 3 : A 3 = u od UQ + u 2 d U2 + • • • + u 2 n d U2n + • • • , 

5 4 : X 4 = U\d ui + u^d U3 + • • • + u 2 rL +id U2ri+1 + • • • , 


S 5 : X 5 


uid UQ + -usd U2 H-+ -— —u 2 n+ id U2 H-, 

3 2n + 1 


S 6 : Xq — uq3 Ui + 3u 2 d U3 + • • • + (2 n + 1) • u 2 n d U2ri+1 + • • • , 


(3.17) 


coo . 
°2n • 


S oo 

2n+l 


x% = d uo „ + 


U 2 n 


(2n + 2)! (2n + 4)! 

(2 n)! “ 2ll+2 2 ! (2 n)l U2n+4 


x : 


2n+l 


= 8n 


+ 


ta„ 


(2ra + 2m)! 

+ m!( 2 n)! U2n+2m '' 

j2n + 5)^ 

I o I / r> , -1 \ I ' ^Uo.t ,J.R I 


U2n+1 ( 2 n + 1 )! “ 2n+3 2 !( 2 n + l)! “ 2n+5 

( 2 n + 2 m + 1 )! 


+ ••• + 


m\( 2 n + 1 )! 


■ t m a 


w 2n+2m.+ l 


+ 


m > 0, 


m > 0, 


where in each case is the infinitesimal generator of the symmetry S/-. These generators form 
a closed Lie algebra as shown in the commutator table below. Since a Lie algebra inherently 
defines a linear vector space, any linear combination of generators (3.17) forms again a sym¬ 
metry generator of the considered system (3.15). And, due to the large variety of symmetries 
admitted, it is possible to combine these symmetry generators (3.17) into several independent 
symmetries such that they are compatible with the posed initial condition given in (3.11). 
Altogether three such compatible symmetries can be constructed from (3.17): 


OO 

Sf : xf = X 2 -J2(k~ 2) 

k =0 

00 

Sf: X* = X 3 + X4-J2 

k =0 

°° / 1 \ 

S 3 : Xl = X 5 + A 6 - £ • u 2k+1 ( 0) • X% + (2k + 1) • u 2 fe (0) • X% +1 ) . 

When generating invariant solutions from (3.18) ^ , all three symmetries independently will 

^The explicit result for the symmetries (3.17) was obtained by augmenting the considered system (3.15) to 
an extended system which also includes all admissible differential consequences of (3.15). Otherwise one runs 
the high risk of performing an inconsistent symmetry analysis. Note that the symmetry corresponding to the 
linear superposition principle has not been included in (3.17), since it cannot be taken to directly construct 
group invariant solutions as we intend to do here. 

*Note that the determination of an invariant function from any of the three symmetries (3.18) is not sufficient 
to be a solution of the underlying system (3.11). In general such an invariant function is endowed with integration 
constants which can only be uniquely determined by plugging the function back into its determining system of 
equations, i.e. here, back into system (3.11). But note that if one of these integration constants can not be 
equated consistently, then the invariant function does not constitute a solution. 


U2k(0) ■ X ^ + 1l2fc+l(0) 


yoo 
^ 2 fc+l 


)■ 


> (3.18) 


U2k(0) 


X%k + “2fc+l(0) • X^fc+l 
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A 2 

A 3 

A 4 

A 5 

a 6 

Voo 

A 2 n 

VOD 

^2n+l 


0 

X\ 

0 

0 

0 

0 

(2n+2)! -\roo 
(2 n)! 2n+2 

(2n+3)! v"oo 
(2n+l)! 2n+3 

X 2 


0 

0 

0 

0 

0 

—(n—2)A^ 

— (n—2)A2^ +1 

X 

CO 



0 

0 

-A 5 

A 6 

-yoo 

A 2 n 

0 

a 4 




0 

As 

-a 6 

0 

VOD 

^2n+l 

A5 





0 

A 3 + A 4 

0 

1 v 00 

2 n+l" X 2 n 

a 6 






0 

-(2n+l)X^ +1 

0 

voo 
A 2 n 







0 

0 

V 00 

^2n+l 








0 


Table 1: Commutator table for the generators (3.17), where [X.;. Xj] = — \Xj,Xj\ := XjXj — XjXi. 


finally give the same unique solution set (3.13), with c n = tt n (0) (3.16), as when solving 
directly the system of equations (3.11). 

4. Example for an infinite forward differential recurrence relation 
— An unclosed system with a non-unique general solution manifold 

In this section we will investigate an infinite system which is not closed; which yet not even 
formally can be regarded as closed. By analyzing the combined system (2.18)-(2.19) for A = 1, 
we are now obtaining a one-dimensional system which, in contrast to the previously in Section 3 
considered system, now features all key properties of the higher level induced (Hopf-equation 
induced) infinite Friedmann-Keller hierarchy of MPC equations (2.13), in particular where the 
crucial property of being a forward recursively defined relation is now included. And it is solely 
this property which now turns the problem into an unclosed problem. 

As it was done in the previous section for the specification A = 0, let us also for A = 1 first 
consider the higher level PDE (2.18) as the following well-posed Cauchy problem 

dtu = d 2 u — x 2 u, for t > 0, with u(x,0) = 4>(x), (4.1) 

where we again assume that the function u = u(x, t ) is decaying sufficiently fast at infinity, i.e. 

lim u(x,t ) = 0, lim d x u(x,t ) = 0, Vt > 0, (4.2) 

x —^iLoo x —>-=boo 

and where also again, for convenience, we want to assume that the initial function cf is nor¬ 
malized to one 

/ OO 

4>(x) dx = 1. (4.3) 

-oo 

To solve this initial value problem (4.1) it is necessary to realize that the following nonlinear 
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point transformation (Polyanin, 2002)1 

t = ^ • (e^ — l) , x = x ■ e 2t , u = u- e~ 2 x2 ~*, ( 4 . 4 ) 

which has the unique inverse transformation 

t = ^ In (1 + 4i), x = - , u = u- ^1 + 4 1 ■ e 2 1 + ii , (4.5) 

4 v ' vTTiT 

maps the original Cauchy problem (4.1) into the following Cauchy problem for the standard 
diffusion equation with constant coefficients: 

dfU = <9| u, for t> 0, with u(x, 0) = (f>(x) ■ e~ 2 x2 . (4.6) 


Important to note here is that the initial time t = 0 as well as the relevant time range t E [0, 00 ) 
both get invariantly mapped to t = 0 and t E [0, 00 ) respectively. Hence, the unique solution 
of the transformed Cauchy problem (4.6) is thus again given by (3.4), but now in the form 


u(x, t) = 


V47rt J-c 






dx ‘, for t > 0 , 


(4.7) 


which then, according to transformation (4.4), leads to the unique solution for the original 
Cauchy problem (4.1) 


u(x,t ) 





(e 2t x-x') 2 

e 4 *- 1 <j>{x') 


-It -' 2 

e 2 


dx', 


for t > 0 . 


(4.8) 


If we again would choose as a initial condition, for example, the Gaussian distribution (3.5) 


(j)(x) 




a, ji E M, a > 0, 


(4.9) 


then the solution (4.8) would have the explicit form 


u(x,t ) = 


7t(ct 2 + e 4t (l + cr 2 ) — 1) 


x 2 e 4t (/x 2 + 2x 2 (1+a 2 )) —4e 2t /IX+/J. 2 

2 2( CT 2 + ,V (1 + ct 2 ) _ 1) ^ t >0. (4.10) 


^This continuous point transformation is not a group transformation, as it neither includes a group parameter 
nor does it include the unique continuously connected identity transformation from which any infinitesimal 
mapping can emanate. 
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Its moments for all t > 0 are 


uniquely given as (shown here again only up to third order) 


n = 0 : uo(t) 



u(x,t ) dx 


n = 1 : u\(t) 



u(x,t ) dx 


n = 2 : v, 2 (t) 



u(x,t ) dx 


n = 3 : us(t) 



u{x , f) dx 


(-l+e 4t ) M 2 


y/2 e 2(l-<T 2 +e 4t (l+CT 2 )) 

\Jl — a 2 + e At (1 + a 2 ) 

(4.11) 

« (-l+e 4t )M 2 

2/x\/2 e 2(i- CT 2 +e 4t (1+<T 2 )) 

\J {1 — cr 2 + e 4t (l + er 2 )) 3 

(4.12) 

f (-l + e 4 t )M 2 

y2e 2(l- CT 2 +e 4i (1 + CT 2)) 


\J {1 — (T 2 + e 4t (l + CT 2 )) 5 


. ( 4e V-(-l + a 2 ) 2 + e 8 *(l + a 2 ) 2 ), 

(4.13) 

n f (-1 +e 4 V 2 

2/i\/2e 2(l-(r2+e 4 *(l + o-2)) 


\J{l - cr 2 + e 4t (l + cr 2 )) 7 


• (4e 4 ^ 2 - 3(— 1 + cr 2 ) 2 + 3e 8t (l + cr 2 ) 2 ). 

(4.14) 


Now, let’s consider the corresponding moment equations when they are independently derived 
from the higher-level Cauchy problem (4.1) (see Appendix D) 

du poo 

—4 = n • (n — 1) • u n _2 — u n+ 2 , with u n (0) = x n ■ 4>(x) dx , n > 0. (4-15) 

tth J -OQ 


Our aim is again to investigate in how far the above infinite system of first order ODEs rep¬ 
resents a well-posed initial value-problem for the moments u n = u n (t). As it was also already 
noticed for the system in the previous section, the first positive observation is that in this case 
too, one can bijectively associate an initial condition u n (0) to each unknown function u n . But 
the next step, deriving the general solution of (4.15), stands in clear contrast to the previous 
section: The system’s general solution is not uniquely specified, i.e. even if sufficient initial 
conditions are imposed, system (4.15) still does not offer a unique solution. In clear contrast, 
of course, to its associated higher level PDE system (4.1), which, as a Cauchy problem, is 
well-posed by providing the unique solution (4.10) with its unique moments (4.11)-(4.14). To 
see this problem explicitly, it is helpful to first recognize that the infinite ODE system (4.15) 
can be rewritten into the following equivalent and already solved formf 


°° l-2i 

U 2 k+2(t) = {~l) k+1 J2 A i ] dt k+l-2i k > 0 ’ 

i =0 

°° jl+2-2j 

w+l(() = (-1)' £ Af w-^=2i 1 a 


(4.16) 


Un the following we agree on the definitions that 


d° _ i di<° 
dip' ■ dtl 


0, and 
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where the coefficients 


A A \k) and A^\l) 


are recursively defined as: 


k—(2i—l) 

Af\k)= ^ ( : 2 k - 2q) ■ (2k - 1 - 2q) ■ A^k - 2 - q ), i > 1, k > 0, 

9=0 

l-(2j- 3 ) 

4 2) (0 = E (2/-2g)-(2/ + l-2g)-^_ 1 (Z-l-g), j > 2, Z > 1, 

9=1 


(4.17) 


with the initial seeds: AqE—1) = 1 , A A \k) = 1, Vfc > 0, and a| 2 ^(0) = 1, A^\l) = 1, Vi > 1, 
respectively. In contrast to the integrated general solution (3.13) of the previously considered 
unrestricted system (3.15), we see that the degree of underdeterminedness in the above derived 
general solution (4.16) is fundamentally different and higher than in (3.13). Instead of integra¬ 
tion constants c n , we now have two integration functions uo(t) and ui(t) which can be chosen 
freely. Their (arbitrary) specification will then determine all other solutions for k > 0 and 
l > 1 according to (4.16). The reason for having two free functions and not infinitely many 
free constants is that the unrestricted system of equations in (4.15) 


= n • (n - 1) • Un- 2 (t) ~ u n+ 2 (t), n> 0, (4.18) 

defines a forward recurrence relation (of order two)^ that needs not to be integrated in order to 
determine its general solution, while system (3.15), in contrast, defines a backward recurrence 
relation (of order two) which needs to be integrated to yield its general solution. 

To explicitly demonstrate that (4.16) is not a unique general solution, we have to impose the 
initial condition u(x, 0) = (j>(x) for the corresponding moments as given in (4.15). The easiest 
and most apparent way to perform this implementation would be in a first choice to choose 
the two arbitrary functions uo(t) and ui(t) as analytical functions, since they straightforwardly 
allow for an expansion into power series 



(4.19) 


where cjf 1 and c ^ are two different infinite sets of constant expansion coefficients. By inserting 
this Ansatz into the general solution (4.16) and imposing the initial conditions as given in (4.15) 
will then uniquely specify these coefficients in a recursive manner as 


Cp ] =0, p < 0; 

4 2) = °> p < °; 


oo 

4 1} = ( _1 ) P • ' U 2p(0) - J2 A i 1] (p ~ X ) • 4-2r> P > °> 

r =1 

oo 

4 2) = ( _1 ) P • “2p+l(0) - E A< S-l(.P) • 4~2r> P > 0- 


(4.20) 


Indeed, the two functions (4.19) with the above determined coefficients (4.20) form the analyt¬ 
ical part of the corresponding unique PDE moment solutions (4.11)-(4.14) when choosing the 
initial-condition function </> as the Gaussian distribution (4.9). In other words, if the unique 
PDE moment solutions (4.11)-(4.14) were Taylor expanded around t = 0, they would exactly 
yield the first four power series solutions (4.19) of the associated infinite ODE system (4.15). 
But, the obvious problem is that the underlying PDE solutions (4.11)-(4.14) do not constitute 


^The order of the recurrence relation is defined relative to the differential operator. 
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analytical functions in the independent variable t, i.e. the Taylor expansions of the unique 
solution functions (4.11)-(4.14) do not converge on the whole temporal domain, but only in the 
following limited range (dependent, of course, on the value of the variance a 2 ) 


0 < <j < 1: 


= 1 : 


<r> 1: 


0 < t < _J 7r 2 + 


0 < t < oo, 


In 


1 — cr 2 
1 + a' z 


0<t<y(2vr) 2 + 


In 


a 2 — l 
a 2 + 1 


(4.21) 


That means, our initial assumption that the first two ODE solutions uo(t ) and ui(t) are 
analytical functions on the global and unlimited scale t E [0,oo), for all values of a 2 , is thus 
not correct. Only for a very limited range, where er 2 ~ 1, this functional choice (4.19) is valid. 
But, if we don’t know the full scale PDE solutions (4.11)-(4.14) beforehand, how then to choose 
these two unknown functions uo(t) and u\(t) for the infinite ODE system (4.16)? The clear 
answer is that there is no way without invoking a prior modelling assumption on the ODE 
system itself. Even if we would choose specific functions fo(t) and f\ (t), which for uo(t) and 
Ui(t) are valid on any larger scale than the limited analytical Ansatz (4.19), we still have the 
problem that this particular solution choice is not unique, because one can always add to this 
choice certain independent functions which give no contributions when evaluated at the initial 
point t = 0. For example, if uo(t) = fo(t) and u\(t) = fi(t), and if both functions /o and fi 
satisfy the given initial conditions at t = 0, then 

2 2 

uo{t) = fo(t) + ■ e~~$, Ui(t) = flit) + i>i(t) ■ e~^, (4.22) 


is also a possible solution choice which satisfies the same initial conditions, where ifo(t) and if \ (t) 
are again (new) arbitrary functions, with the only restriction that, at the initial point t = 0, 
they have to increase slower than e 7 o/* and e 77 ?* respectively. 

Note that uo(t) and u\(t) are not privileged in the sense that only these functions can be 
chosen arbitrarily. Any two functions u n *(t) and u m *(t ) in the hierarchy (4.18) can be chosen 
freely, where n = n* and n = m* are some arbitrary but fixed orders in this hierarchy such 
that n* / m* + 2k, V&: E Z. For example, for the choice n* = 2 and m* = 3 in (4.18) its 
underdetermined general solution will take the form 


uo{t) = u 0 ( 0 ) - [ u 2 {t')dt', 
Jo 


OO 

2=0 

Ui(t) = 7x1(0) - [ u 3 (t')dt\ 
Jo 


/ ft \ d k ~ 2i 

do,k+i-2i Uo(0) - ^ u 2 {t')dt'\ - dtk _ 2i u 2 (t) 


, k > 1 , 




J = 1 


rt \ d^^— 2 ^ 

5o,l+2-2j (7Xi(0) - J ^ u 3 {t')dt'j - dtl+1 ^ 2j U 3 (t) 


l > 2 , 


(4.23) 

where now, compared to the alternative general solution (4.16), the functions u 2 {t) and u 3 {t) 
are the unclosed terms, instead of uo(t) and ui(t). 

That the infinite system (4.18) is underdetermined and thus unclosed is also supported when 
performing an invariance analysis on it. To transform system (4.18) invariantly, complete 
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arbitrariness exists in that three arbitrary functions are available in order to perform the 
transformation: One for the independent variable t and two for any arbitrary but fixed chosen 
dependent variables u n * and u m *, such that n* / m* + 2k, Vk E Z, i.e. where in effect the 
transformation of two functions u n * = u n * (t) and u m * = u rn * ( t ) can be chosen absolutely freely. 
Note that the outcome of such an (unclosed) invariance analysis does not result into symmetry 
transformations, but only into weaker equivalence transformations which invariantly only map 
between unclosed systems (see also the discussion partly done in the Introduction; for more 
details on this issue, see Frewer et al. (2014a, b) and the references therein). 

Hence, when choosing, for example, n* = 0 and m* = 1, then infinitely many functionally 
independent equivalence transformations of the following form are admitted by system (4.18) 


: X = id t + r] UQ d U0 + ij ui d ui 

+ Vui (£, Vuo , r) ui , u 2 , u 3 ) d U2 
+ Vu 3 (£, riuo , Vui , u 2, u 3 ) d U3 

+ ••• + Tj U2n ^,T] uo ,T] ul ,U2,U 3 , . . . ,U2m'^2n+l)9 U 2 n 

+ Vu 2n+1 {£,,Vu 0 ,r] ul ,U2,U 3 , . . . ,U2n,U2n+l)d U2rl+1 + ■■■ , (4.24) 


where £ = £(i,«o,ui), r] Uo (t,uo,ui) and r) Ul (t,uo,ui) are arbitrary, free choosable functions. 
Consequently, the infinite set of invariant transformations (4.24) do not form a closed Lie 
algebra (in contrast to the infinite set (3.17) of the previously considered system). In par¬ 
ticular, when imposing the initial condition u(x, 0) = <j>(x) for the corresponding moments 
as given in (4.15), then infinitely many and functionally independent invariant (equivalence) 
transformations can be constructed which all are compatible with this arbitrary but specifi¬ 
cally chosen initial condition. Because, since e.g. the infinitesimals £ = £(t,uo,ui), r] Uo (t, uo,u\) 
and r] Ul (t,uo,ui) can be chosen arbitrarily, one only has to guarantee that each of the initial 
conditions u n (t = 0) = u n ( 0), for all n > 0, gets invariantly mapped into itself. This is achieved 
by demanding all infinitesimals to satisfy the restrictions 


Vn 0 u 0 , Ui. 


£(t,u 0 ,u i) 
= 0 


{t=0;u n =u n (0) ,V)i>0} 


= 0 , 


{t=0;u n =Mn(0),Vn>0} 

^]u2 n (C) ) Vui ) ^2 j ^3 j ■ ■ ■ j ^2n j ^2n+l) 

Vu2n+1 (£ ! Vuo J Vui , U2 , u 3 , . . . , U 2 n i u 2n+ 1) 


r] Ul (t,u 0 ,ui) 

=n n (0),Vn>0} 

{t=0;u n 


= 0, n > 1 

{t=0-,u n =u n (0),Vn>0} 



= 0, n > 1 

{t=0;Un=n,i(0),Vn>0} 



= 0 , 


(4.25) 


(4.26) 


where only the three infinitesimals £, r) UQ and r/ ui can be chosen freely, while the remaining 
infinitesimals 2 n and V 2 n+i, f° r all n > 1, are predetermined differential functions of their 
indicated arguments. The conditions (4.25), in accordance with (4.26), can be easily fulfilled 
e.g. by restricting the three arbitrary functions £, r] uo and r] ui to 


£(£, uq, ui) = f 0 (t, u 0 ,ui) ■ e , 


Vu o (t,u 0 ,ui) = go(t,u Q ,ui) ■ e (-o-“o(o)) 2 , r] Ul (t,u 0 ,ui) = ho(t,uo,ui) ■ e (-i— 1 (°)) , 


(4.27) 


where /o, go and ho are again (new) arbitrary functions, however, now restricted to the class 
of functions which are increasing slower than e 1 ^ at r = 0, where 


r = yj t 2 /tf + (uq - u O ( 0 )) 2 / 7 | + («i - ^i( 0)) 2 /7 l- 


( 4 . 28 ) 
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And, since in this case all differential functions r] U2n and rj U2n+1 , Vn > 1 have the special 
non-shifted affine property 


•iu 2n 


{?=0; T]u 0 =0; r] ui = 


(4.29) 


the conditions (4.26) all are automatically satisfied by the above restriction (4.27). Hence, 
an infinite set of functionally independent (non-privileged) invariant solutions for each mo¬ 
ment u n = u n {t), Vn > 0 can be constructed from (4.27), which all satisfy the given initial 
condition u n (t = 0) = u n ( 0) as imposed in (4.15). 

That a priori no unique general solution, nor any unique invariant solution can be con¬ 
structed, provides the reason that the PDE induced ODE system (4.18), although infinite in 
dimension, has to be treated as an unclosed system. It involves more unknown functions than 
there are determining equations, although formally, in a bijective manner, to each function 
within the hierarchy a corresponding equation can be mapped to. But, since the hierarchy 

(4.18) can be equivalently rewritten into the solved form (4.16) (or into (4.23), etc.) it explic¬ 
itly reveals the fact that exactly two functions in this hierarchy, e.g. uo(t) and u\{t), remain 
unknown, and without the precise knowledge of their global functional structure all remain¬ 
ing solutions u n fit ) for n > 2 then remain unknown too. Even when posing sufficient initial 
conditions is not enough to yield a unique solution for the lower-level ODE system (4.15) as 
it is for the higher-level PDE equation (4.1). Without a prior modelling assumption on the 
ODE system (4.15), this system remains unclosed. Fortunately, the solutions of this particular 
case (4.11)-(4.14) possessed an analytical part in their functions for which the assumed Ansatz 

(4.19) expressed the correct functional behavior, though only in a very narrow and limited 
range. But, of course, for more general cases such a partial analytical structure is not always 
necessarily provided, and an Ansatz as (4.19) would then be misleading.! 

Hence, the infinite forward recursively defined system of order two (4.18) does not possess 
a unique general solution. The degree of arbitrariness in having two unknown functions cannot 
be reduced, even when a sufficient number of initial conditions are imposed, simply due to 
the existing modus operandi (4.22) when constructing possible valid solutions. The same non¬ 
uniqueness problem we face for invariant solutions due to the construction principle (4.27), in 
that infinitely many functionally independent equivalence transformations can be constructed 
which all can be made compatible to any specifically chosen initial condition. Hence, the 
existence of an invariant solution for such a system is without any value, since it just represents 
an arbitrary solution among infinitely many other, equally privileged invariant solutions. 


5. Conclusion 

At the example of a lower-dimensional ODE-system this study has shown that an infinite 
and forward recursive hierarchy of differential equations carries all features of an unclosed 
system, and that, conclusively, all admitted invariance transformations may only be identified 
as equivalence transformations. To obtain from such systems an invariant solution which should 
show a certain particular functional structure is ultimately without value, since infinitely many 
functionally different and non-privileged invariant solutions can be constructed. In order to 
obtain valuable results, the infinite system needs to be closed by posing modelling assumptions 
which have to reflect the structure of the underlying (higher-level) equations from which the 
infinite system emerges. 

^Note that not only an analytical Ansatz as (4.19) would be misleading, but also any other, functionally 
different Ansatz can be misleading too, if no further external (exogenous) information is at hand. For example, 
choosing instead of an analytical a Fourier-expansion Ansatz would globally fail as the underlying functions 
(4.11)-(4.14) do not show any time-periodic behavior, at all. 
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Although in this study only an example of a lower-dimensional system has been considered, 
it is well to be expected that the very same non-uniqueness issues also occur in the higher¬ 
dimensional Friedmann-Keller hierarchy of moments for turbulent flows, and that it would be 
instructive to remember this example when trying to generate solutions or invariant solutions 
from the infinite and unmodelled Friedmann-Keller equations as its motivated e.g. in the work of 
Oberlack et al. Besides the problem of being defined forward recursively, the infinite Friedmann- 
Keller hierarchy faces another and further substantial problem, namely the singular problem of 
not showing a symmetric and thus complete coupling among all equations. As it was discussed 
at the end of Section 2, the problem that in each equation the occurring unknown higher order 
moment is not directly coupled to the corresponding next higher order equation, additionally 
even increases the degree of underdeterminedness of the Friedmann-Keller equations, which 
thus undoubtedly renders it altogether to an unclosed system. 

Hence, the approach and the conclusions given in Oberlack & Rosteck (2010); Oberlack & 
Zieleniewicz (2013); Avsarkisov et al. (2014); Waclawczyk et al. (2014); Oberlack et al. (2014), 
about deriving physically relevant invariant solutions as so-called turbulent scaling-laws from 
first principles for the lower level moments within the infinite Friedmann-Keller hierarchy, is 
heavily misleading; all the more so as these first principle invariant solutions are additionally 
based on physically inconsistent symmetries (for more details, see Frewer et al. (2014a, b)). 
The general argument that these turbulent scaling laws are ’’clearly validated” in Oberlack & 
Rosteck (2011a,6); Oberlack & Zieleniewicz (2013); Avsarkisov et al. (2014); Oberlack et al. 
(2014), by comparing to direct numerical simulations (DNS), is based on a fallacy (see Section 5 
in Frewer et al. (2014a)). The problem is that this ’’validation” is always only performed for the 
lowest order moments, which, of course, can be robustly matched to the DNS data since there 
are enough free group parameters available to be fitted; in the very same sense as John von 
Neumann famously remarked many years ago (as quoted by Enrico Fermi to Freeman Dyson, 
see e.g. Dyson (2004)): 

'’With four parameters I can fit an elephant, and with five I can make him wiggle his trunk.” 

But, as soon as the higher order correlations functions get fitted, the curve-fitting procedure of 
Oberlack et al. consistently fails. This outcome is supported by a mathematical proof given in 
Frewer et al. (2014a) (see Appendix D), which clearly shows that the Lie-group based turbulent 
scaling laws for all higher order velocity correlations as derived in Oberlack & Rosteck (2010); 
Oberlack & Zieleniewicz (2013); Avsarkisov et al. (2014); Waclawczyk et al. (2014); Oberlack 
et al. (2014) are not consistent to the scaling of the lowest order correlation function, which for 
shear flows is the mean velocity field. In other words, the proof in Frewer et al. (2014a) shows 
that for the lowest correlation order n = 1 no contradiction exists, only as from n > 2 onwards 
the contradiction starts, i.e. while the mean velocity field can be robustly matched to the DNS 
data, it consistently fails for all higher order correlation functions and gets more pronounced 
the higher the correlation order n is. 

Hence, by always only matching or comparing the DNS data to the lowest order velocity 
correlation of a specific flow, as it was done particulary in Oberlack & Rosteck (201 la,6); Ober¬ 
lack & Zieleniewicz (2013); Avsarkisov et al. (2014); Oberlack et al. (2014), i.e. by matching 
or comparing only to the mean velocity profile for wall-bounded (shear) flows and only to the 
two-point velocity correlation for isotropic flows, is definitely not enough to be a true validation 
of the Lie-group based scaling theory. In particular, as this theory by Oberlack et al. which 
specifically considers the infinite (unclosed) system of all multi-point correlation equations and 
which thus is designed and laid-out to be a scaling theory for all higher order velocity cor¬ 
relations, special attention has to be devoted to the prediction value of all those correlation 
functions which go beyond the lowest order. And exactly this has been investigated in Frewer 
et al. (2014a) (Section 5), which then gives a completely different picture than the ’’validation” 
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procedure in Oberlack & Rosteck (2011a,6); Oberlack & Zieleniewicz (2013); Avsarkisov et al. 
(2014); Oberlack et al. (2014) is trying to suggest. 

Taken together, the key reason for this failure is twofold: i) Oberlack et al. considers the 
infinite and forward recursive Friedmann-Keller hierarchy of moments as a closed system, which, 
as was clarified in this study at an analogous lower-dimensional system of ODEs, is not the 
case; ii) Oberlack et al. considers all layers of statistical description of turbulence endogenously 
without involving at the same time the underlying deterministic Navier-Stokes system — a 
methodological approach which obviously is incomplete and even incorrect, because it are the 
deterministic equations which due to their spatially nonlocal and temporally chaotic behavior 
induce the statistical equations and not vice versa. In other words, in order to obtain valuable 
and physically consistent results within a statistical description of the physically quantifiable 
world, it is necessary to include all underlying processes up to the fluctuating level into this 
description and not to ignore them if they exist. 


A. Alternative form of a transverse vector field 


The transverse component A J_ (x) of an arbitrary vector field A(x), when decomposed into a 
longitudinal and transverse component A(x) = A 11 (x) + A J ~(x) respectively, and, which decays 
sufficiently fast for |x| —>• oo, is defined as (see e.g. Stewart (2008)) 

i , . . . f V' • A(x') 

A ± (x) := A(x) + V 


TT^x', 


(A.l) 


47t|x — x' 

confirming its defining (solenoidal) property V • A J_ (x) = 0, when making use of the standard 

1 


identity 


<5 3 (x-x') = -A- 


47t|x — x'| 

Relation (A.l) can be written in an alternative form, by rewriting it as 


(A.2) 


A x (x) = A(x) + V 


= A(x) + V / V 


V ' ■ A (*') 

47t|x — x'| 

, ( A(x') 


\47t|x — x' 

= A(x) — V / A(x') • (V 


1 


= A(x) + V / A(x') 


V 


= A(x) + / A(x') ■ V® V 


47t|x — x' 

1 

4-7r|x — x'| 

1 


47t| x — x'| 


d 3 x' - V J A(x') • ( V' 
d 3 x 
d 3 x' 

dV 


47t|x — x'| 


d 3 x' 


= / d 3 x' d 3 (x - x') • 1 + V ® V- 


1 


A(x'), 


(A-3) 


47r|x — x'| 

where we made use of Gauss’ theorem for the second term in the second line, which resulted 
into a vanishing surface integral due vanishing contributions of the vector field A(x) at infinity. 
Hence, the alternative relation (A.3) for a transverse vector field allows to define an operator 
T 7 of orthoprojection to solenoidal vector fields: A e-?• V[A] = A ± . The explicit expression 
for "P acting on any arbitrary vector-field A vanishing sufficiently fast at infinity then has the 
form given by (A.3) (see e.g. also Appendix B in Frewer et al. (2014a)): 


V[A\= / d 3 x'p(x - x') • A(x'), 


(A.4) 
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with the kernel 

p(x - x') := <5 3 (x - x') • 1 + V <g> V—- 

y y ’ 47t|x-x' 

By construction, the projection properties of 'P are 

V[A] =A, if V • A = 0, 

V[A] =0, if V x A = 0, e.g. if A = V<p. 

B. The Hopf equation in two different representations 

If we assume that V- y / 0, for any arbitrary vector field y = y(x) which is decaying sufficiently 
fast at infinity, the functional Hopf equation as given by (2.1) can be represented in two alter¬ 
native ways. Either in a representation relative to the full vector field y, or in a representation 
relative to its projected transverse (solenoidal) part y^. 


(A.5) 

(A.6) 


B.l. The Hopf equation in the full-field representation 
According to the defining relation (2.2), the functional Hopf equation (2.1) will take the form 



d 3 x 


d 


5 2 4 > 


x i 


+ vA- 


(5<E> 


dxp Sy a (x)6yp{x) 8y a fx) 


cTx + i / <p(x 


d 2 


< 5 2 4 > 


dxadx /3 5y a {x.)Syi3(x.) 


d 3 x 


x i 


. d 


5 2 <& 


a \ ,3 

+ z^A-—— ) d 3 x 


+ i 


dxp 6y a (x)5y/3(x) 8y a (x.) 

d 1 \ d 2 


5 2 $ 


dxly 47t|x — x'| / dx a dxp 5y a (x.)5y/3(x.) 


d 3 x'd 3 x 


d 


5 2 <$> 


x i 


A \ ,1 

+ vA -—— | ePx 


dxy 8y a (x)8yp(x) 8y a (* 


+ i 


' 'll fxl 

( 9 1 ) 

/ Vay*-) 

\dx a 47r x — x' ) 


d 2 


5 2 $ 


d 3 x'd 3 x. 


(B.l) 


Note that in going from the second to the third equality, the incompressibility condition in 
functional space (Hopf, 1952, p. 97)^ has been made used of 


d 8 4* 

dx a 5y a (x.) 


(B.2) 


^Note that the incompressibility condition only refers to the velocity field u of the flow, and not to the 
auxiliary field y. In other words, although we assume V ■ y ^ 0, i.e. y ^ y x , we still have incompressibility in 
the velocity field V ■ u = 0, i.e. u = u x ; and in functional space it is expressed as (B.2). 
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B.2. The Hopf equation in the transverse representation 

In order to transform the Hopf equation (2.1) into a representation where only the transverse 
fields appear, it is necessary to first transform the functional derivatives according to relation 
( 2 . 2 ) in applying the corresponding functional chain rule 


<fya(x) 

<5 2 


= / (T z 


3 „ <5 2/k(z) 6 


5y a (x) 5y^(z) ’ 


(B.3) 


5y a (x)Sypf 


= td\ 


■3„ j3„/ S Vk( z ) Sy\( z ') 


X 


<fya(x)<5y,s(x) <5y, 


-L ( 


+ / d 6 zd 6 z 


5y a (x) %(x) 5y 3 r(z)5y^(z') 


rv ( B - 4 ) 


Let’s first consider the transformation of the single functional derivative in the viscous term: 


<%«(*) 


= J d z l 5 aK 5 6 ( z - x) + — 


d f d'<S 3 (x' -x)^,^ 5 


47t|z — x'| 


<%K X ) 

5 

^Va (x) 
£ 

<%K x ) 


- / d 3 z 


- / d 3 Z 


+ / d 3 z 


_ 9 _ 

dz K 

d d 


5 3 (x'-x; 


a 


^( z ) 


d 3 x' 


<9x' Q 4 tt|z - x'| J 5y^(z) 


dz K dx a 47t|z - x| / 5y^-(z) 


5 


dx a 47t|z - x| / <9z K 5y^-(z) 


5 


(B.5) 


If this operator now acts on the characteristic functional 4*, which itself transforms invariantly 
(Hopf, 1952, p. 93) 

$[y(x),i] = ^[y L (x),t], (B.6) 

then the second right-hand-side term in the last line of (B.5) results to zero due to the incom¬ 
pressibility condition (B.2). Hence we obtain the following invariant result 


(54>[y(x),t] _ < 5$[y- L (x),i] 


6y a {yd) 

Similar for the second derivative in the inertial term, which also gives the invariant result^ 


(B.7) 


5 2 4> 


5y a (x)$yp(- 


X 


= J d 3 z crz I S aK 5 6 (z - x) + — 


5 /^(x'-x)^ 


5p X 5 3 (z , -x) + - t 


47t|z — x'| 
8 r <^'<5 3 (x"-x) 


= / d 3 zd 3 z 5 nK 5 3 (z — x) — 


dz x J 47r|z'— x /y 
d d 1 


d 3 x" 


5 2 $ 


SyH z )Sy x (z') 


cLc a 47 t|z — x| 


5p X S (z - x) - — 


9 5 


1 


J 2 $ 


94 dxp Att\z’ - x| / 5y^(z)<5y 3 -(z') 


^Note that since the functional relation (2.2) for y ± is linear in y, the first term of the chain rule (B.4) gives 
no contribution. 
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5 2 $ 


( x ) 

+ [ d 3 zd 3 z 5 aK 8 3 (z-x)- i , ± 

J dxp 47t|z' — x| 5y£(z) \dz' x 5yf{z , ) / 

+ J d 3 z d 3 z d/3\5 3 (z r — x) 

d 


d 


d 


1 


d 5$ 


d 6$ 


+ / d 3 z d 3 z! 


dx a 47r|z — x| (z') \&Zk <%/«(*) / 

1 \ f d 1 \ <9 5 / <9 


9x Q 47r|z — x| I l dxp 47r|z' — x| J dz K Sy^( z) l 5j/^(z') 


%a(x)^(x)’ 

where again we made use of the incompressibility condition (B.2). Now, since the characteristic 
functional $ transforms invariantly (B.6), the transformed Hopf equation (2.1) in the transverse 
representation finally takes the form 


<9$ 

~dt 



(. d 5 2 <S> M> \ 3 

\dxpSy^(x)5yj(x) +U %H X ) J 


(B.9) 


C. The MPC equations in two different representations 

Similar as in the case for the Hopf equation (see Appendix B), the MPC equations for the 
incompressible Navier-Stokes velocity field u = u(x, t) can be represented in two alternative 
ways too. Again, either in a representation relative to the full MPC fields^ 

Hct { n } := H ai an := (u ai (xi) ■ ■ ■ u an (x n )), (C.l) 

or in a representation relative to the transversely projected MPC fields 

H a {n} ■= : = («ai(xi) ■ • •U a „(x n )) ± . (C.2) 

Note that due to (2.10), along with its evaluation (2.11), we have the obvious invariant relation 
= This invariance can be readily verified independently from relation ( 2 . 10 ), by 

recognizing that i) the (fluctuating) velocity held is solenoidal d ak u ak (xk ) = 0 , for any index k, 
i.e. u at (xfc) = u^ k (x*.) , and ii) that any usual (non-functional) differential operator acting 
locally on any point x& is commuting with the averaging (ensemble) operator (•). Hence we 
have the relation 


da k (Uai(xi) ■ ■ ■ U ak (xk) ■ ■ 


= (u Q1 (xi) • • 

• d ak u ak (x fc ) ■■■u oln (x n )> = 0 

(C.3) 

(u Ql (x i) • 

" ««„(**)} 

= (u ai (xi) • 

" u ci„( x n)) ■ 

(C.4) 


But careful: When taking the smooth and regular limit of a zero-correlation length between 
any two points |x*, — x/| —>• 0, the invariance (C.4) no longer holds, because obviously 

d ak lim (u ai (xi) ■ ■ ■ U ak (*k) ■ ■ ■ U ai (*l) ■ ■ ■ Ua n (-Xn)) 

= (u ai (xi) • • • Ua k (xk) ■ ■ ■ da k U ai (xfc) • • • U an (x n )) ^ 0. (C.5) 


Mhe //-notation for the (instantaneous) MPC functions in (C.l) is taken from Oberlack & Rosteck (2010); 
see also Frewer et al. (2014a) for a critical discussion on using this notation, in particular for invariance analysis. 
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C.l. The MFC equations in the transverse representation 
The underlying MPC-generating equation is given by (2.12) 


94> n 

~df 


Va( x ) 


+ i/A- 


d<W 


d <5 2 4> n+1 
dx P dy^{x)8yj{x) ' " * ^Va ( x ) 


d 3 x, Vn > 1, 


(C. 6 ) 


where the hierarchy of functionals 4> ra is given by (2.9), along with (2.10) and (2.11). Three 
terms are contributing to equation (C. 6 ): the temporal term on the left-hand side and the 
inertial and viscous terms on the right-hand side, which for each order n evaluates to (shown 
here only up third order) 


n = 1 : 


n = 2 : 


n = 3 : 


^f = f d 3 xiy«i(xi) (4(^ 1 (xi)) ± ) , 


<5 2 $ 2 


2 • 1 


5y x {x)5yf(x.) 2 ! 

Ah 1 


(u a (x)ug(x)) , 


<^a( x ) 

a $ 2 

“dT 


= * («q(x)) , 


= J d 3 x 1 d 3 x 2 y^ 1 {x 1 )y^ 2 {x 2 ) ( _ ^?^ ti «i( Xl )“«a( x 2))' L ) , 


5 2 4 > 3 


%a( x )^( x ) 

<f $ 2 


/<*3 V 


3-2 


<%K x ) 

9<h 3 
“dT 


= J d 3 x 2 y^ 2 (x 2 ) (-|j (u a (x)« Q2 (x 2 )) ± ^ 


= J d 3 x i d 3 x 2 d 2 x 3 (xi) y^ 2 (x 2 ) y„ 3 (x 3 ) 



(C.7) 


(C. 8 ) 


(C.9) 

-« 2 ( x 2 ))^, 

(C.10) 

)u a3 (x 3 ))^ , 

(C.ll) 


(C.12) 
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_ *3! 9t^“ 1 ( Xl ) W "2( X 2)«a3( x 3))" L ) , 


<5 2 $ 4 


%a( x )^ 2 /y( x ) 


= J d 3 x 3 d 3 x 4 y^ 3 (x 3 )y^ 4 (x 4 ) 


4-3 j_ ■ 

—jj- (u Q {x.)up (x) ri Q3 (x 3 ) ^Q!4 ( x 4)) ) , 


d<l> 3 

^ 2 /q ( x ) 


(C.13) 


(C.14) 


= J d 3 x 2 d 3 x 3 y^ 2 (x 2 )y^ 3 (x 3 ) ^-i ^ (u a (x)u a 2 (x 2 )u a 3 (x 3 )) ± ^ , (C.15) 


and so on for all higher orders. When inserting these results for each order into equation (C. 6 ), 
two things should be noted or to be carried out before dropping the integrals in order to obtain 
the MPC equations in the transverse representation: i) The integrals for the inertial and viscous 
terms must be fully symmetrized in their integration variables in accord with the integrals from 
the temporal term, e.g. consider the second order result for the inertial term (C.ll), which, 
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through equation (C. 6 ) and in accord with result (C.10), has to be symmetrized to 


J d 3 x d 3 X 3 i/a (x) ?/4 (x 3 ) ^ <u Q (x)^(x)u a3 (x 3 ))- L 

= j d 3 xi d 3 x 2 2/4(xi) 2/4( x 2 ) ^ ^ <u ai (xi)u Q2 (x 2 )ri fTl (xi))- 


2 0 x 2^2 

= [ d 3 xi d 3 x 2 2/4 ( x i) 2/4 ( x 2) I * ~ ^ H, 


1 <9 

+ X --(u Q , 1 (xi)u Q2 (x 2 )u t72 (x 2 )) 


2 9xi jC 


OL\OL 2 (J\ 


i d fT 

+ - - H, 


2 dx 


2 ,<72 


OL\OL2^1 


, (C.16) 


and ii) the invariance condition tp = H n (C.4) is only valid if the MPC function of n-th order 
= (#4.. an ) is also evaluated at n distinct points x n . As soon as it’s evaluated in less than 
n distinct points, the invariance no longer holds (C.5). These lower-dimensional MPC functions 
will be denoted by H,4 as it was already done in (C.16), in order to avoid any uncertainties in 
the notation. For them, as already said, the invariance condition does not hold, i.e. H,^ ^ H n . 
Following both these instructions i) and ii), we finally obtain the hierarchy of MPC equations 
in the transverse representation (shown here only up third order) 


n = 1 : 

n = 2 : 

n = 3 : 


dt 


ox i iCri 


d 2 H aia 


U 11 Ct\Ot2 _j_ U11 OL\OL2(Tl _|_ U11 OL\OL2(T2 ^ ^ n -OL\OL2 ^ — j.— 

dt dxi t(Tl dx 2 ,a 2 dxi^dxi^ dx 2 , U 2 dx ; 

Q TT TJ- L TT _L tX-L 

( ^' Jrl aia2Q , 3 | u 11 (\\ 02CV,3(T1 _|_ u 11 ot.ia.20L3(T2 ___QlQ2Q3^3 

dt dx i >cri dx 2)(T2 dx 3)a 3 

, d 2 H a j , 
dx 3)l73 c 


rfi TJ U 

^ C/ ^^Q;ia2Q!3 ^ ^ ^ z aia2Q:3 ^ 

9xi j(T 1 3a:i i(T1 dx 2 )(T 2 dx 2 ^ 2 


(C.17) 

(C.18) 


(C.19) 


where the coupling within the hierarchy is given through the definition 


TJ -L 

11 a.\.. .OLncr k 


lim H, 

a+l ^X/e! ^ra+1 


Q!i...a n Q; n +i 


Vfc between 1 < fc < n. 


(C.20) 


Note that the infinite hierarchy of MPC equations is accompanied by the following obvious 
incompressibility constraints 


dH ai ... ari 

dxk,a k 


\/k between 1 < k < n, 


£)TJ- L 


<9^, 


= 0 , 


«fe 


Vfc,/ between 1 < (k,l) < n , for k ^ l. 


(C.21) 

(C.22) 


C.2. The MPC equations in the full-field representation 

The full-field representation of the MPC equations follows directly from (C.17)-(C.19) by 
decomposing each of the inertial transverse terms according to relation (2.2). But, as writ¬ 
ten, relation (2.2) only applies for vector fields. Therefore it is necessary to first extend this 
definition to also include tensor fields of arbitrary rank such that it applies to (C.17)-(C.19). 
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Following the principle of (2.2) it can straightforwardly be generalized to construct transverse 
(solenoidal) tensor fields 


A 


at 


,(x) = A 


ai—ot n 


(x) + 


n 


d', A n 


n 


k= 1 




47r|x — x'| 


ft 3 x', 


(C.23) 


where A n = (A ai ... an ) is an arbitrary tensor field of rank n decaying sufficiently fast at infinity. 
By construction, (C.23) is in accord with (2.2) since its full-ranked divergence vanishes^ 

= 0. (C.24) 


It is obvious that relation (C.23) also further extends to mufti-point tensor fields, hence we 
obtain the following decompositions of the transverse MPC fields appearing in (C.17)-(C.19) 
(shown here only up to second order) 


d#ai CT i( x i) = dH aiul (xi) + 3 r d 3 x[ (x() 

dxi )(T1 3xi, CT1 dxi , ai J 4tt|xi -x(| 3x' q , 3x' )(7 , 

= ggWxi) + 3 r ^(^KK'K)) 

dxi :<Tl dxi , Ql 7 47 t|xi-x / 1 | 

_ 3# q , 10 - 1 (xi) 3(p(xi)) dH aiai (y.j) g%](xi) 

3xi j(Ti 3xi, Qi ' dxi )CTl dxi )OL1 


n = 2 : 


^q 1 q 2CT1 ( x ^ x 2) fl#q iaag2 (xi,X 2 ) 

3xi lfri 3^2, (72 

3g aia2CT1 (xi,x 2 ) ! 1 3 r d 3 x.[ 3 2 ff a / Q2fT /(x(,x 2 ) 

3xi, CT1 2 3 xi , ai J 47t|x 1 -x' 1 | dx\ ^ dx\ 

1 3 r d?yf 2 3 2 ^ aia ^/(xi,x^) 

2 3x 2 ,o 2 J 47t|x 2 -x' 2 | ^x i 2q ,Qx' 2(j ^ 

+ 3# aia2CT2 (xi,x 2 ) + 1 3 r ft 3 x(, 3 2 ff aia ^ a ./(xi,X 2 ) 

3x 2 , 0-2 2 3x 2 ,a 2 j 47r|x 2 -X2| ^ X 2 ,a' 2 ^ X 2 ,a' 2 

1 3 r ft 3 x( 3 2 ff a / iQ2 , J /(x / 1 ,x 2 ) 

+ 2 3xi, Ql J 4tt|xi - x'J 3x' l Q , 3x( 

__ 3ff Qia20 - 1 (xi,x 2 ) 3(p(xi)u a2 (x 2 )) dH Ql g 2a2 (xi,x 2 ) 3(u ai (xi)p(x 2 )) 

3xi )(Tl 3 xi, Q1 3x 2i0 - 2 3 x 2 ,o 2 

3iftp 1 Q 20 - 1 (x X , X 2 ) 3/q 2 [i](xi, X 2 ) 3i7q 1 q 2tr2 (xi,X2) 3J ai [2]( x l, x 2) 

3xi, CT1 3xi, ai 3x2,02 3x2,02 

(C.26) 


^Note that only a full ranked divergence on A x 1 ... oln (C.23) vanishes, i.e. when all indices are contracted. A 
partial divergence where not all indices are contracted does not result to zero. Hence, the partial divergences of 
the transverse MPC fields appearing in (C.17)-(C.19) are non-zero. 
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wich then according to (C.17)- 
in the full-field representation 

d 2 H ai 


n = 1 : 


re = 2 : 


n = 3 


5/ [i] 

-—- V - 

9xi, CT1 9xi, CT1 


= 0, 
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?presentation (C.17)-(C.22), the above equations 


— 1,U1 “ ™ -L,C»1 

As before in the transverse field representation 
through the definition 


0, 


(C.28) 


(C.29) 
are coupled 


H ai ...a n a k = lim H ai ... an a n+ 1; VA: between 1 < k < re, 

Xn+l^Xfc, OL n -\-\—>&]<. 


and are accompanied by the incompressibility constraints 


dH ai ... an _ ^ between 1 < k < n, 


9x 


k,a k 


dlrvi ••• 


Ql 101+1 ——M = o, \/k, l between 1 < ( k , l ) < n, for k ^ l. 

<JXk,a k 


(C.30) 


(C.31) 

(C.32) 


Note that the above system agrees with the MPC hierarchy of equations derived in Oberlack 
&; Rosteck (2010), from which also the above notation has been borrowed (see also Frewer 
et al. (2014a, b)). These are i) the instantaneous (equal-time) multi-point velocity correlation 
functions of order n > 1 


H a{n} = H ai ... an (x i,... ,x n ,f) = (u Ql (x!,t) ■■■u otn (x n ,f)), (C.33) 

where H ai ... anCTk (xi,..., x n , x*,, t) is a lower dimensional moment which emerges from the next 
higher one H ai ... an0lri+1 (xi,..., x n , x ra+ i, t) in the limit of zero correlation length |x n+ i—x*.| —>• 0 
along with a corresponding change in the tensor index a n +i —>• <7*,, for all k within the range 
1 < k < n (C.30), and then ii) the instantaneous (equal-time) multi-point pressure velocity 
correlation functions of order (re — 1) 


4 { „_ 1} [J] = («a 1 (x 1 ,t)---Ua I _ 1 (xj_ 1 ,t) •p(x i ,t)-u ai+1 (xi +1 ,t)---u an (x n ,t)), (C.34) 

where p is the instantaneous pressure field, which in the unbounded domain evolves according 
to Poisson’s equation as (see e.g. McComb (1990)) 



d£pu a {x.',t)up(x',t) 
47r|x — x'| 


(C.35) 
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D. Derivation of the ODE moment equations from the underlying PDE 

In a similar fashion as the infinite (lower level) Friedmann-Keller hierarchy of moments (2.13) 
emerge from the underlying (higher level) Hopf equation (2.6), the analogous lower level ODE 
moment equations (2.19) are induced by the higher level PDE (2.18) in the following way: 
Multiplying the unbounded PDE for u = u(x, t) 

dtu = d%u — A • x 2 u, (D.l) 

by x n , for all n > 0, and integrating over M by assuming that partial integration is justified, 
i.e. by assuming the natural boundary conditions 


lim u(x,t) = 0, lim d x u(x,t) = 0, Vt > 0, 

x —>-zboo x —>-=boo 


(D.2) 


then one obtains for each order n the following infinite system of coupled ODEs 

= n ■ (n - 1) • u n - 2 — A • u n+ 2 , Vn > 0, 


dt 


where the moments u n = u n (t) are defined by 

/ OO 

x n ■ ufx,t) dx, Vn > 0. 

-oo 

Formulating (D.l) as a Cauchy problem, by including the initial condition 

/ OO 

f>(x) dx = 1, 

-OO 


(D.3) 


(D.4) 


(D.5) 


where <f> is some arbitrary but integrable function, it will correspondingly restrict the infinite 
system of equations (D.3) to a first order ODE initial value problem, in that it has to satisfy 
the restrictions 

POO 

u n { 0) = x n ■ 4>(x)dx, Vn > 0, with uq(0) = 1. (D-6) 
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